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Abstract 

In this paper, an argon-like canonical system is studied. We introduce 
five hypothesis to deal with the total potential of the system. Then the 
balanced liquid-gas coexistence phenomenon is analyzed. Good equations 
of state and phase diagram are given. 
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1 INTRODUCTION 

Much attention has been paid in recent years to the hard core Yukawa (HCY) 
potential as a model for the pair interactions of fluids [1]. The hquid state theo- 
ries such as the Mean Spherical Approximation (MSA) [2] and the Self Consistent 
Ornstein-Zernike Approximation (SCOZA) are proposed. Recent studies of the 
HCY fluid can be found in [2,3] and references therein. Worthvifhile similar 
hard-core Sutherland potential 

+oo,nj < 2ro; 
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was less investigated, i and j sign two particles, r^j is their distance, and tq 
is interpreted here to be the radius of a hard core. 5 is a constant. When 
^ = 1, this potential is Sutherland potential. It is generally accepted that ^ > 1. 
In this paper, wo mention out five hypothesis to deal with the total potential 
energy of a balanced canonical system, on the basis of which we obtain the 
canonical partition function by integrating under the help of Eq.(l). 

In section 2 of this paper, the equation of state and the chemical potential 
are gotten by thermodynamical formulas. So do their reduced forms. Phase 
diagrams under five similar hard-core Sutherland potentials are illustrated . 

In section 3 of this paper. We conclude that there is no proper similar hard- 
core Sutherland potential which can bring completely right interpretation to 
the balanced liquid-gas coexistence phenomenon by the theory in this paper. 
We find that a similar Sutherland potential with its exponent ^ = 1.8 leads to 
more accurate forecast not only to the critical coefficient of argon but also to 
the phase diagram and the relation of pressure versus temperature. 

2 THEORY 

2.1 The Canonical Partition Function 



The case of a balanced argon-like gases canonical system in 3 dimension space 
is considered firstly. The average volume v of each particle is defined as 

where n is the particle number density, V is the volume of the system, and N 
is the particle number of the system. 

Suppose that volume t; is a three dimensional sphere with its radius vd- 
Then we have 

3 '^^ - iV - n- 

It is natural that ro > tq- They are illustrated clearly in Fig.(l). 

As far as the potential energy is considered, we mention out such hypotheses 

as 

1. The potential energy uij of two particles with their distance rij is described 
byEq.(l). 

2. U = NU^. 

U is the total potential energy of the system and is divided to each particle 
equably. Ui is called the total potential energy of particle i, which benefits from 
other (A'^ — 1) particles and is the sum of the potential energy contributions of 

[N — 1) particles to it. 

Q rr. . _ Tj. . — Itt. . 
— *j — — — 2 "^^ ' 

Uj^i is the potential energy contribution of particle j to the total potential 
energy of particle i, and Ui^j is the potential energy contribution of particle i 
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to the total potential energy of particle j. The potential energy Uij is divided 
to these two particles equably. Thus, 

Ui = {Ui^i + U2^i + ... + U^i-i)^i + U(^i+i)^i + C/(i+2)^i + ... + t/jv^i) 

= I (Uu + U2i + ... + + i7(i+i)i + ?7(,+2), + ... + Um) . (4) 

4. When the thermal fluctuation is omitted: Ui^i + U2^i + ... + + 

When the idea that the total potential energy is divided to each particle 
equably is accepted, the total potential energy of arbitrary particle i keeps the 
same in the case we discussed. Our case is a balanced argon-like gases canonical 
system. When the thermal fluctuation is omitted, we do think that the sum of 
the potential energy contributions of (iV — 1) particles to particle i is equal to 
the potential energy contribution of particle j to particle i in the universe space 
of this balanced system, ri is the down limit of the integration and r2 is the up 
limit of the integration. 

5. ri = Td, T2 = +00. 

The down limit is supposed to be r^. And the up limit is supposed to be 
+CX). From Eq. (1), we know that ri should take the value of ro, which is 
the radius of the hard core, ri = ro is equal to such an idea that integration 
Ii = J^" UijAnr^jdrij <^ I2 = J^^ UijAnrfjdrij . Hence Ii is not taken into 
account in the calculation of potential, rzj > ro is illustrated in Fig.(l). 

From hypotheses (1-5), we have the result below 

U = NUi=N {Ui^i + U2^^ + ... + U^i-l)^i + U(^i+l)^^ + ... + UN^i) 

= Nlj^Uijd%. (5) 

We define 

Uij = / Uijd^rlj, (6) 
Jv 

and integrate it as follows 

Uij = / Uijd^r\j 
Jv 

r+00 

= / Uij^irrjjdrij 
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3-3^ 

Here we applied r/j > ro, ^ > 1. 
The Hamiltonian H reads 



TV 



2m 

where m is the mass of one particle, pi is the momentum of particle i. 
Thus canonical partition function can be solved as 



1 



exp{-pH)d'''yd'''\, (9) 



with (3 = l/ksT, h Plank constant. We calculate Eq.(9) as follows 
Q = T77TW / exp{-pU)d'^r, 



1 



mx 



3N 



zn{T,V) (10) 



with thermal wavelength A = -- — ? is the Boltzmann constant and 

Zi\r{T,V) is called position partition function normally. 
Inputting Eq.(5) to Eq.(lO), we get 

^i^/exp(-,^ii^(..)3-).3.. 
= ^ / exp(/?B'7Vn(-i+«)d3^r 



^ exp(/3B'7Vn")V/^, (11) 



3JV 



with a = ^—1, B' = — ^r^(^)'^~^' and Vj is the free volume, n is t/ie particle 
number density which keeps conservation when a balanced canonical system is 
considered. And we do not consider the fluctuation of n in this paper. For the 
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existence of the radius tq of the hard core, each particle excludes a volume 
and hence for N particles the excluded volume is ^^^iV. The free volume, 
therefor Vf = V - Nb with b = Thus Eq. (11) reads 

Q = — ^ eM-PB'Nn") {V - Nbf . (12) 



2.2 The Equation of State 

We can get the equation of state as follows 

'SlnQ 



P = kBT 



N,T 



dV 

aB'rf+^. (13) 



V -Nb 

When we choose cr = 1, the equation of state is written as 

(P + B'n^) {v-b) = UbT. (14) 

This is just the form of the VDW equation of state. 

Here we sign the pressure of the gases Pi, the pressure of the liquids 
the chemical potential of the gases yu^, the chemical potential of the liquids 
^2, the critical temperature T^, the critical pressure Pc, the critical particle 
number density ricthe particle number density of the gases ni, the particle 
number density of the liquids n2, the reduced temperature T*, the reduced 
particle number density of the gases n* , the reduced particle number density of 
the liquids nj, the reduced pressure of the gases P^* ,the reduced pressure of 
the liquids Pj*, the particle number of the gases A^i, the particle number of the 
liquids A^2 , the volume of the gases V\, and the volume of the liquids V2. The 
relations between them are 

T* = T/Tc, (15) 
ni = 7Vi/yi, (16) 
n2 = N2/V2, (17) 
n* = ni/ric, (18) 
nl = n2/nc, (19) 
Pi* = Pi/Pe, (20) 
P2* = P2/PC. (21) 
At critical point, the function P = P{N, V, T) has such qualities as 

|^k=0, (22) 
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^k=0. (23) 
Thus we get the critical data by solving Eq.(22) and Eq.(23). They are 

/cBre = a(a + l)(^)2s'<, (25) 

Pc = :^B'nl+\ (26) 



nckeTc ^ 4(g + 1) 
Pc (t{(t + 2) ■ 



When cr = 1, they read 



(27) 



ne = ^, (28) 

M = (29) 

^c = ^, (30) 

C = ^^2^ = 8/3. (31) 

We are very familiar with these results which are just the reduced data from 
the VDW equation of state. C is the critical coefficient. When a = 0.7432, 
C = 3.4201. 

Inputting Eq. (24-26) to Eq.(13), we get the reduced equation of state 

^ An*T* (g+l) _ (a + 2) n*("+^) 

{{a + 2)-n*a)a a ' ^ ' 

2.3 The Chemical Potential 



Here we can solve the chemical potential as follows 

'ainQ' 



M = -ksT 



ON 



V,T 



= -B' (. + 1) + fc.T m j + ^ - ^ m (33) 

Inputting Eq. (24-26) to Eq.(33), we get the critical chemical potential. Then 
the reduced chemical potential is calculated to be 

M = -, , (34) 
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where Ai , A2 and A3 are given by 



Ai 



In- 



n a 



AaT* 



+ 



A2 = 



b{a + 2-n*a)J {a + 2y {a + 2 - n*a) {a + 2) 



6TV StttoB'TV^'^+i) (cr + 1) 



(a + 2r 



6- (a + 2)"+' /i2 



, /, (7 a 3, 87rmB'cr('^+i)(cr + 1)\ 4ct 

-1 + In \ In 77 

\ 26 2 2 b-{a + 2f+^ {a + 2f 



(35) 
(36) 

(37) 



2.4 Balanced Liquid-gas Coexistence Canonical System 



We consider the case of the balanced liquid-gas coexistence canonical argon-like 
system in 3 dimensional space, which is a unit-two phases system. Here we do 
think the two phases are described by the same canonical partition function 
as Eq.(12). There are three balanced conditions which must be satisfied when 
these two phases are balanced. They are thermal condition 

Ti = T2, (38) 

dynamic condition 

Pi = P2, (39) 

and phase condition 

/xi = /X2. (40) 
For Eq. (15-21), these three conditions can be expressed in the reduced unit 

as 

n = n, (41) 

P* = P2, (42) 
= M2- (43) 

Prom Eq.(32), we have 

^ 4^Tr(a + l) _ (a + 2)nt^-+^) 
^ {{a + 2)-nla)a (7 ' ^ ' 

^ 4n:,Ti (g + l) _ (^ + 2).*'-+^^ 
^ {{a + 2)-n*2a)a a ' ^ ' 



Prom Eq.(34), we have 
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M2 = 



(^l)2 + (^2)2 



(^3)2 



The solution of Eq.(41) and Eq.(42) is 

(K)(-+^) - (1 - n,K)6) (1 - n^jnm 

' ' (^ + l)(42)'(K)-("2)) 



(47) 



(48) 



The solution of Eq.(41) and Eq.(43) is 



In 



("l)(l-"c(";)fa) 



-("l)fc 



{n2)(l — nc(n^)b) (l — ric (n]^ )6) "^T— nT^n^^bJ 

Thus we get a function /(n|, 712) = easily from Eq.(48) and Eq.(49). 

((nf)(-+i) - (n^)(-+i)) (1 - n,{nm (1 - ne(n^)6) 



(49) 



/K."2) = 



(a+l)((nl)-(n^)) 
((nt)- - (n^)-) 



Jjj ("l)(l-"<:("2)^) ^ «c(«^)b ric(n|)6_ 

(n|)(l-rac(nt)b) (l-nc(n^)&) (l-nc(n|)b) 



(50) 



If (T is chosen to be a constant, can be gotten in the way of numerical 
computation by computer when an arbitrary n* is fixed. Then T* is obtained 
easily from Eq.(48) or Eq.(49). And P* is solved by Eq.(32). Table.(l) is the 
theoretic data when a = 0.8. Fig. (2) is the correlation of n* and under 
different similar hard-core Sutherland potentials signed by different a. Fig. (3) 
is the phase diagrams corresponding with different a. Fig. (4) is the curve of 
InP* versus 1/T*. In Ref.[4], we have introduced one method called polynomial 
approximation to simulate the relation of two variables. Here we can get the 
relation of InP* versus 1/T* by this method when a is fixed. In Ref.[5], the 
theoretic relation of InP* versus 1/T* was given by this method in the case of 

(7=1. 

In 1945, E.A.Guggenheim collected the data of the balanced liquid-gas co- 
existence system from experiments and gave out the correlation of pi, P2 and 
T* by the empirical equations[5,6] below 



For that 



we have 



p\ = l + 0.75(1 - T*) - 1.75(1 - T*)^/^ 
p*2 = l + 0.75(1 - T*) + 1.75(1 - T*Y'^. 

p\ = mni/mric = ni/ric = n^, 
P2 = mn2/mnc = n2/nc = n^, 

= 1+ 0.75(1 - T*) - 1.75(1 - T*)i/^ 



(51) 
(52) 

(53) 
(54) 

(55) 
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n* = 1 + 0.75(1 - T*) + 1.75(1 - T*^/^. (56) 

In Eq.(51-52), is the reduced density of the gases and P2 is the reduced 
density of the Hquids. As far as argon system is concerned, the inaccuracy 
of these two equations is generally only one or two parts per thousand of P2 
or of pI when T* > 0.60Tc[5,6]. So, it is acceptable to consider the data of 
{T*,nl{T*),n2{T*)) from Eq.(55-56) as the experimental ones in this tempera- 
ture region[5]. 

E.A.Guggenheim gave a numerical analytic result of the relation between 
the reduced temperature and the reduced pressure from experiments by equa- 
tion[5,7] 

P* = exp(5.29 - 5.31/T*), (57) 

which best fits the experimental data for argon when T* > 0.60Tc except a tiny 
region near the critical point [5, 7]. Thus it is acceptable to consider the data of 
{P*,T*) from Eq.(57) as the experimental ones, too[5]. The data from Eq.(55- 
57) are illustrated as experimental data in Fig. (2-4) and Table. (2) to compare 
with the theoretic ones deduced above. 



3 CONCLUSIONS AND DISCUSSIONS 

In this paper we advanced five hypothesis to deal with the total potential energy 
of a balanced canonical system. Hypothesis 1-4 are easy to be accepted. The 
key is hypothesis 5. Actually ri = Td, ^2 = -|-oo is a simple approximation to a 
concrete sytem. In labs, the volume of the sytem we consider is finite:r2 7^ -|-oo. 
Additionally, only when we consider a static sytem in 3 dimensional space, can 
the average volume of a particle be regarded as a sphere. And rzj is effective, 
ri = rn is only a simple effective approximation. 

On the basis of the five hypothesis, we got the canonical partition function 
Eq.(12). Thus all the thermodynamic quantities will be solved by Eq.(12). 
Following, we analyzed the balanced liquid-gas coexistence canonical argon- 
like system. The results are partly illustrated in figures and tables. Fig. (3-4) 
indicates that a in the region (0,1) leads to more perfect forecast to experimental 
data. Applying numerical calculation by computer, we find a = 0.8 is the best 
among the five values of a chosen in this paper only when phase diagram and 
the relation of In P* versus 1 /T* are considered together. It is clearly illustrated 
in Figure. (3-4). Then the similar hard-core Sutherland potential is fixed to be 
the form of 



+00, nj <ro, 



In this case, the equation of state is 

Now the critical coefficient C is equal to 3.2143. But a = 0.8 is not the best 
when the relation of nl versus n2 is considered. Figure. (2) suggests that a = 0.9 
and 1 are better. 
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Now wc will ask whether a proper cr, which can bring completely right 
forecast to experiments when the liquid-gas coexistence phenomenon is con- 
sidered, exisits or not by the theory in this paper. Suppose that such a a = 
(To exists. Thus Eq.(50) will be right when arbitrary two experimental data 
terms (n|,n2) are considered. From Eq. (55-56), we get two experimental data 
terms:(0.1963, 2.0137) |t*=o.8600, (0.3928, 1.6822) |t»=o.950o- We input them to 
Eq.(50) and get the numerical value of a by Matlab software. The result is 
o'|t*=o.8600 = 0.98 72 ^ cr|T*=o.950o = 1.0883. Thus we conclude such a proper 
a does not exist. But it does not mean that there exists no value for a corre- 
sponding with the potential explaining the experimental data properly. A new 
theory may be do. Actually, when the liquid-gas phase transiton to second order 
is considered, our work does not work well for the obvious fluctuation near to 
the critical point, which is ommittcd in hypothesis (4). We will discuss it in de- 
tails in future work. Of course, we see that this work offered the proper critical 
coefiicient C = 3.4201 of argon by a simple equation of state with a = 0.7432. 
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Figure 1: Two radius 



Figure 2: The curve of n| versus 712- The real one is from Eq.(60-61). The 
rest are all from Eq.(55) with different a: square-cr = 2; diamand — a = 1.5; 
star a = 1; dot — cr = 0.9; cross-tr = 0.8. The star one is the result of the 
VDW equation of state. 



11 



Figure 3: The phase diagrams. The real one is from Eq.(60-61). The rest are 
from the theoretic data with different a: square-a = 2; diamand — cr = 1.5; 
star — a = 1; dot — a = 0.9; cross-cr = 0.8. The star one is the result of the 
VDW equation of state. 
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Figure 4: The curve of InP* versus 1/T*. The real one is from Eq.(62). The 
rest are theoretic data with different a: square-cr = 2; diamand — cr ~ 1.5; 
star — a = 1; dot — a = 0.9; cross-c = 0.8. The star one is the result of the 
VDW equation of state. 
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nt 0.0106 0.0292 0.0535 0.0851 0.1266 0.1827 0.2627 0.3928 

2.8097 2.6574 2.5328 2.4122 2.2860 2.1458 1.9785 1.7521 

P* 0.0158 0.0483 0.0931 0.1521 0.2287 0.3280 0.4590 0.6411 

T* 0.4911 0.5743 0.6360 0.6906 0.7426 0.7944 0.8484 0.9085 

Table.l : Theoretic data when a = 0.8. 

nl 0.0106 0.0292 0.0535 0.0851 0.1266 0.1827 0.2627 0.3928 

2.5894 2.4958 2.3965 2.2899 2.1734 2.0423 1.8873 1.6822 

P* 0.0284 0.0562 0.1007 0.1670 0.2599 0.3840 0.5434 0.7412 

T* 0.6000 0.6500 0.7000 0.7500 0.8000 0.8500 0.9000 0.9500 

Table.2 : The data from Eq.(55-57). 



14 




The reduced density of the liquids 
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